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ABSTRACT 


A molecular beam time-of-flight technique is studied as 
a means of determining surface stay times for physical 
adsorption. The experimental approach consists of pulsing a 
molecular beam, allowing the pulse to strike an adsorbing 
surface and detecting the molecular pulse after it has 
subsequently desorbed. The technique is also found to be 
useful for general studies of adsorption under non- 
equilibrium conditions including the study of adsorbate- 
adsorbate interactions. 

The shape of the detected pulse is analyzed in detail 
for a first-order desorption process. For mean stay times, 
r, less than the mean molecular transit times involved, the 
peak of the detected pulse is delayed by an amount approxi- 
mately equal to r. For r much greater than these transit 
times, the detected pulse should decay as exp(-t/r). How- 
ever, for stay times of the order of the transit times, both 
the molecular speed distributions and the incident pulse 
duration time must be taken into account. 

Estimates of r were obtained from the experimental 
results for Xe, Kr and C0 2 on nickel and for Xe on copper 
surfaces as a function of the surface temperature T g . Xe 
and C0 2 were found to have mean stay times of about 10”^ to 

10‘ 3 sec over a range of T from 125 to 105°K, Observed 

s 
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values of t for Kr were about 10 J to 10 sec in the range 

99 to 92°K. The effect of He, Ar and Xe ion bombardment was 

to reduce the stay time for Xe on nickel by a factor of 10 

to 100 at a given temperature. This is attributed to burial 

of neutralized ions near the surface. 

Binding energies, E 0 , and pre exponential factors, r , 

are estimated by fitting the results to the expression 

T * r exp(E /kT). Values of r are found to be about 10”*^ 
o o o 

sec in order-of-magnitude agreement with theoretically pre- 
dicted values for a localized adsorption model. Values for 
E 0 /k ranged from about 2600°K for Kr on nickel to about 
3250°K for Xe on nickel. 
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CHAPTER 1 


INTRODUCTION 


1.1 Adsorption and Stay Times 

The adsorption of molecules by solid surfaces has been 
the subject of a large number of investigations for many 
years. Molecular adsorption from the gas or liquid phase 
plays an important role in technical processes such as the 
purification of gases and liquids, dyeing of materials, and 
surface catalysis. Adsorption from the gas phase, in par- 
ticular, is an essential mechanism for the operation of 
cryogenic vacuum pumps. However, adsorption is not always 
a desirable process. The presence of even a monolayer of 
adsorbed gas can sometimes lead to significant degradation 
of a device whose performance is strongly dependent on its 
surface properties, e.g., thermionic emitters and optical 
coatings. Designers of vehicles operating in the space 
environment as well as designers of ultra-high vacuum 
systems are especially concerned with preventing contamina- 
tion of critical surfaces. Thus, a strong interest exists 
in the adsorptive behavior of a variety of gas-surface 
combinations under conditions where the adsorbate coverage 
constitutes one monolayer or less. 

Adsorption may be defined as that process which causes 
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the molecular density near the surface to be higher than 
the density in the gas phase. Adsorption is classified 
according to the nature of the forces which cause it. If 
the adsorption is caused by weak, van der Waals forces, it 
is known as physical adsorption or, simply, physi sorption . 

If the adsorption involves the transfer or sharing of elec- 
trons between the gas and solid atoms, it is then known as 
chemical adsorption or, simply, chemisorption . Physical 
adsorption is characterized by the fact that the resulting 
binding energies between the gas and solid are usually of 
the same order of magnitude as the heat of liquefaction of 
the gas. 

Adsorption, in general, is the result of competition 
between the adsorption rate (rate of deposition) and 
desorption rate (rate of evaporation). When the adsorption 
rate exceeds the desorption rate, the surface coverage 
increases with timet conversely, when the desorption rate is 
greater, the surface coverage decreases. Of course, once 
the gas-solid system has achieved a state of thermodynamic 
equilibrium, the adsorption and desorption rates must be 
equal and the surface coverage remains constant. 

The majority of experimental studies of adsorption have 
been carried out under conditions where the gas-solid system 
was allowed to reach a state of equilibrium. In most cases, 
the approach to equilibrium was quite rapid and adsorption 
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was considered to occur spontaneously^^. While this type 
of study gives useful information on equilibrium surface 
coverage, more information can be obtained about the adsorp- 
tion history of a given molecule from non-equilibrium 
experiments. For example, the time and surface coverage 
dependence of adsorption and desorption probabilities can be 
determined, in principle at least, from a detailed knowledge 
of the adsorption and desorption rates. 

One approach to obtaining the non-equilibrium informa- 
tion is to measure the mean stay time (sometimes called 
dwell time, residence time, or sticking time) that adsorbed 
molecules spend on the surface. The idea that a molecule 

striking a surface spends a finite time on the surface was 

( 2 ) 

first proposed by Langmuir' and later was developed by 
(*1 \ 

Frenkel w/ from statistical mechanics. The stay time con- 
cept has been examined at length by deBoer' ' to show its 
usefulness in understanding both the equilibrium properties 
of adsorption and the mechanisms which govern the rate at 
which equilibrium is reached. Therefore, knowledge of the 
stay time is of considerable value in the study of adsorp- 
tion. 


1 , 2 Stay-Time Measurement 


In order to determine the requirements for measuring 
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the stay time, it is necessary to examine the expected 
behavior. For a first-order desorption process (desorption 
flux proportional to surface coverage), the Frenkel model 
predicts that the mean stay time is proportional to 
exp(E./kT_) where E n is the binding energy, k is the 
Boltzmann constant, and T s is the temperature of the solid 
at its surface. For physi sorption, where E Q is comparable 
to the heat of vaporization, stay times become measurable 
only as the temperature approaches the liquefaction temper- 
ature of the gas which, for many gases of interest, is below 
room temperature. Also, based on the predicted exponential 
dependence, the stay time will increase rapidly with 
decreasing temperature. While large stay times (of the 
order of seconds) are more easily measured, they in turn 
lead to large surface coverage. Thus, if low surface 
coverage is of interest, stay-time magnitudes must be 
limited, 

1.3 The Present Experiments 

The purpose of the present experiments was to examine 
a method suited to the measurement of 9tay times greater 
than a few microseconds. The technique was one in which a 
thermal energy, molecular beam was pulsed and the pulse of 
molecules was directed onto the surface. The subsequent 
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desorption rate was detected and the time evolution of the 
detected pulse was used to infer a mean stay time. For this 
technique, molecular transit times were important, but they 
were measured as part of the experiment and their effects on 
the shape of the detected pulse were taken into account. 

The experiments were conducted over a range of surface 
temperatures such that the stay time varied by at least an 
order of magnitude about the mean molecular transit time. 
Methods for estimating the stay time from the detected 
pulse were evaluated over this range of stay times. The 
effects of varying the gas species, the incident-pulse 
duration, the incident molecular flux, the surface material 
and means of preparation, and ion bombardment of the sur- 
face have been studied. 



CHAPTER 2 


THEORY AND BACKGROUND 


2.1 Theory 

( 2 ) 

As mentioned in the introduction, it was Langmuir' 
who first proposed that a gas molecule striking a surface 
would spend a finite time before departing. Prior to his 
proposal, it was thought that a critical temperature T 

V 

existed such that, for T>T , molecules were reflected 

c 

instantaneously and, for T < T , atoms were permanently 

C 

adsorbed. 

The stay-time concept was first developed quantita- 
( 1 ) 

tively by Frenkel'-" using the methods of classical 
statistical mechanics. He assumed thatt 

(1) an individual gas molecule is bound to a 
surface with energy E Q r 

(2) the number of molecules adsorbed is suffi- 
ciently small that interactions among 
themselves are negligible » 

(3) adsorbed molecules retain two translational 
degrees of freedom parallel to the surface t 

(4) internal degrees of freedom can be ignored i 

(5) an adsorbed molecule is bound as a simple 

one -dimensional harmonic oscillator vibrating 
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perpendicular to the surface with frequency 
u 0 (period t q = 27 t/u) <3 )\ 

(6) the gas-solid system is in thermodynamic 
equilibrium. 

Frenkel’s derived result for the mean stay time t isi 

T = t q exp(E 0 /kT) (2.1.1) 

where T is the temperature of the gas-solid system and k is 
the Boltzmann constant. This relation has been widely used 
to predict mean stay times as well as to infer values of t q 
and E q from experimental observations of the temperature 
dependence of r. 

The physical significance of the mean stay time may be 
illustrated by considering the desorption of an ensemble of 
molecules by a first-order rate process. If a is the num- 
ber of molecules adsorbed per unit area at any time t and 
P d is the probability of desorption per unit time, then 

d o _ p . 

a — -Jr , o 

TF d a (2.1.2) 

If P d is constant, Eq. (2.1.2) may be integrated to give 


<r a (t) = * a (0)exp(-P d t) 


(2.1.3) 
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The mean stay time for the ensemble is 


r = 



dt 


(2.1.4) 


Substituting for a (t) from Eq. (2.1,3) and integrating 

a 

yields 


t = i/P d (2.1,5) 

This interpretation of the mean stay time is usually valid 

as long as the surface coverage is much less than a mono- 

layer. At higher surface coverages, adsorbate-adsorbate 

interactions become significant and the process is no 

longer first order. In this case, a mean stay time is still 

defined by Eq. (2.1.4), but it is meaningful only when the 

coverage history is also specified. 

Although the Frenkel relation, Eq. (2.1.1), is widely 

used, more general predictions of theoretical stay times 

(4 o 

have been presented by deBoer' * The initial development 
of these predictions is similar to that of Frenkel and 
follows directly from consideration of steady-state 
conditions, i.e., where the desorption rate equals the 
adsorption rate. If F^ is the incident flux and P & is the 
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adsorption probability (sometimes called sticking prob- 
ability), it follows for first-order desorption that 

% ■ Vi T (2 - 1 - 6) 

The flux F^ is usually known from the state of the system 
(e.g., temperature and pressure). For example, for a gas in 
equilibrium with the surface, the kinetic theory of gases 
predicts 

F i = n e(2iSr) (2.1.7) 

for a eas of number density n and molecular mass m , 

£ & 

Substituting this into Eq. (2.1.6) and solving for r, one 
obtains 


T 



( 2 . 1 . 8 ) 


For a system of volume V and adsorbing surface of area A , 

s 



(2.1.9) 
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where N and N are the total number of molecules in the 
a g 

adsorbed phase and gas phase, respectively. Eq. (2.1.8) 
can then be written as 



This relation provides a convenient starting point for the 
prediction of t by the methods of equilibrium statistical 
mechanics. 

The evaluation of r from Eq. (2.1.10) has been treated 
(3) 

by deBoer w using classical statistical mechanics. The 
results for mobile and localized adsorption are of parti- 
cular interest. Mobile adsorption refers to the case where 
molecules are bound to the surface but move freely in direc- 
tions parallel to the surface. Localized adsorption refers 
to the case where molecules are bound so strongly in local 
sites that their translational motion parallel to the sur- 
face is negligible. The development for these two cases is 
given in Appendix I, and a simple extension of deBoer’s 
approach is proposed which includes some quantum-mechanical 
effects. 
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For a mobile adsorbed gas, the result is 


T 


e j/4t 

slnh(VW T <= eXP< V M > 


( 2 . 1 . 11 ) 


where $-q is the Debye temperature of the solid. In the 
high-temperature limit (T»0 D ), Frenkel's result is 
obtained. However, for T<0p, the predicted values for t 
are substantially less than those predicted by the Frenkel 
expression. 

The result for localized adsorption is given by 

47rh 2 <7 m (0JJ/2T) 2 

■g M D [zsinMfljjADp T o e *P<V kT) (2-1>12) 

where is the monolayer surface coverage and 2jrfi is 
Planck's constant. The high-temperature limit for this 
expression is of course altogether different from the 
Frenkel expression. 

The quantitative differences in stay times predicted 
by Eqs. (2.1.11) and (2.1.12) may be illustrated by writing 
these equations in the form 


r » exp (E^kT) 


(2.1.13) 
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Thus, for a given binding energy and temperature, one need 

only investigate the behavior of the preexponential factor 

r' , Sample calculations given in Appendix I for xenon on 

nickel at 100°K show that t' q for mobile adsorption is only 

-1 3 

slightly less than r p both being of the order of 10 ■' sec. 

However, for localized adsorption, t q is of the order of 
IQ" 17 sec. Therefore, stay times are expected to depend 
considerably on the type of adsorption process. 

The general form of Eq. (2.1.10) allows a natural 
extension of the theory to include such effects as restric- 
tion of internal and rotational degrees of freedom as well 
as a more exact quantum-mechanical treatment of bound 
states . 

An analysis of first-order non-equilibrium adsorption 
has been carried out by Pagni^*^ using a classical mole- 
cular collision model and non-equilibrium statistical 
mechanics. His results for steady-state adsorption point 

f 

out the explicit dependence of r Q on the binding energy. 

While the time-dependent case is also formulated, no solu- 
tions for t are presented. 

2.2 Previous Experiments 

A number of stay-time measurements have been made 
prior to the present work. Clausing^ 11 ^ , in the first 
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attempt to measure stay times, directed a stream of Cd 
vapor onto a rotating disk maintained at a constant temper- 
ature. The atoms then desorbed and were collected by a 
condensation plate. It was expected that the mean stay time 
could be inferred from the angular displacement of the con- 
densate. However, the stay times encountered in these 

experiments were less than the time-resolution of the appar- 

-6 

atus (-10" sec). Clausing later devised a technique 
where the time for a gas pulse to diffuse through a capil- 
lary tube^ 2 ^ was measured. Stay times as large as ~10”^ 
sec were measured for A.r, Ne and N ? on glass. Recent exper- 
iments of this same general type have been conducted for He 

( * o ) 

on Cu near liquid helium temperature' J . 

Recent measurements of stay times have mainly involved 
the use of modulated molecular beams^^* 22 ^ . Most of these 
studies' ~ ' have been conducted for chemisorbed gases. 

For all such measurements, a molecular beam is interrupted in 
some manner by a shutter prior to striking the surface of 
interest. If the shutter is opened (or closed) in a time 
that is short compared to the stay time, the desorption rate 
is that for a step-like increase (or decrease) in the depo- 
sition rate. If the beam is pulsed (shutter opened and 
closed) and if the pulse time-width is small compared to the 
stay time, then the desorption rate is that for an instan- 
taneously depose + ed ensemble. For either type of fast 
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shutter behavior, the desorption rate for a first-order 
process is simply a function of the mean stay time. Other- 
wise, the desorption rate depends on the deposition history. 
Also, if shutter-to-surface transit times are not negligible, 
they must be taken into account. The desorbing molecules 
are usually detected by some sort of flux- or density- 
sensitive detector placed in the desorption path. If 
surface-detector transit times are not small compared to 
times of interest, then these too must be taken into account. 
In only two of these studies 20 ^ were shutter duration 
and transit-time effects evaluated in detail. In the other 
cases, they were either negligible or were accounted for 
approximately. 

Stay times have also been measured by the technique 
commonly known as flash desorption. The technique and its 
results have been reviewed by Ehrlich^ 23 ). Although it is 
used primarily for determining binding energies, it has 
been used to estimate chemisorption stay times. The disad- 
vantage of this technique is that the temperature dependence 
of the stay time must be known a priori. 

Physisorption stay times have not been studied nearly 
so thoroughly as those for chemisorption. This is due in 
large part to the experimental difficulties. The present 
study investigates these difficulties and, under selected 
conditions, resolves them to a point where meaningful data 
are obtained. 



CHAPTER 3 

EXPERIMENTAL APPARATUS 

In "this chapter, the molecular— beam and vacuum system, 
the time-of-f light detection system, and the target and 
target preparation techniques are described. The arrange- 
ment of the various components is shown schematically in 
Fig. 3-1.1. The geometry of the beam-defining elements and 
the time-of-f light (TOF) detection system is shown in Fig. 
3.1.?. The vacuum system and instrumentation and the stay- 
time apparatus are further described by the photographs of 
Figs. 3.1.3 and 3.1.4. 

3.1 Molecular Beam and Vacuum System 

The experiments were conducted in a stainless-steel 
vacuum chamber evacuated by a liquid-nitrogen-trapped, oil- 

diffusion pump. The chamber pressure was typically 

-8 

2 x 10" torr with no molecular-beam flow and about 
2 x 10' 7 torr with the beam on. (All test chamber pressures 
reported herein were measured using an ionization gauge 

calibrated for nitrogen, and the quoted values are as-read 
gauge values. ) 

High-purity (nominally 99 . 995 %) gas was inlet to the 
molecular beam source chamber through a variable leak valve. 
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Figure 3.1.1 Schematic of the Experimental Apparatus 
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DETECTOR POSITION A 
( All dimensions in centimeters ) 



*1 * 0.32 


DETECTOR POSITION B 

Figure 3.1.2 Geometry of Beam-Target-Detector for 

Incident Beam and Stay-Time Measurements 
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Figure 3.1*3 Photograph of the Vacuum Chamber and Instrumentation 
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Figure 3.1.^ Photograph of Stay-Time Apparatus 
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The source chamber was maintained at ambient temperature 

(« 295°K) and the source -chamber pressure p Q was measured 

with a capacitance-type electronic manometer. 

The molecular beam was formed by a multichannel source, 

a photograph of which is shown in Pig. 3.1.5* This source 

is composed of a large number (wl.7 x 10^) of quartz 

capillary tubes fused into an array and in turn fused in the 

center of the end plate of a 1.9 cm diameter glass tube. 

Each individual channel has a nominal diameter of 5 x 10“^cm 

and a length of 0.15 cm. The channels are assembled in a 

hexagonal pattern with a nominal overall diameter of 0.1 cm. 

The molecular beam thus formed is then collimated by an 

orifice as shown in Pig. 3.1.2. 

The characteristics of the molecular beam formed by 

this source were the subject of a separate investigation^ 

and are described in some detail in Appendix II. Prom 

measurements of the beam centerline flux and speed obtained 

in that study (with the detector in position A of Fig. 

3.1.2), the flux and energy of the beam striking the target 

(as shown by position B of Fig. 3.1.2) can be estimated. 

For example, for Xe gas with p Q * 10 torr, the intensity of 

the beam at the target is estimated to be 3.4 x 10*^ 

-2 -1 

atoms • cm - sec . The measured mean energy per atom for 
Ar was about 0.064 eV (some 25 % higher than that for an 
effusive beam at the source temperature); Xe is expected to 
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Figure 3.1.5 Photograph of Multichannel Source 
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have the same mean energy as Ar since Ref. (24) indicates 
that, for the source pressures used in the present experi- 
ments, the multichannel source is operating in the limit of 
isentropic expansion. Although the angular distribution of 
beam flux was not measured, based on Ref. (?5), it is esti- 
mated that 75 % of the molecules diverged less than +6° from 
the beam axis. For completely free molecular flow through 
the beam defining elements, the maximum beam divergence 
should not exceed +4°. 

3.2 Time-of-Flight Detection System 

The time-of-flight (TOP) detection system was designed 
after that of Hagena, Scott and Varma^ 2 ^ and is shown 
schematically in Fig. 3.2.1. The system includes a chopper, 
a detector, and detector signal processing instrumentation. 
The chopper is a thin 15.2 cm diameter disk with 0.1 6 cm 
radial slots cut 0,64 cm deep at 90° intervals around its 
periphery and is driven by a synchronous motor. For the 
stay-time experiments, the chopper rotational speed is held 
constant at about 1200 rpm. Direct beam speed distributions 
were determined at a chopper speed of about 3000 rpm to 
obtain reasonable time-of-flight resolution ^ 2 ^ . As a slot 
passes across the molecular beam, a molecular pulse passes 
through and the "start time" for this pulse is sensed by a 
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Figure 3.2.1 Schematic Diagram of the Time-of-Flight 
Instrumentation 
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light-photodiode combination located di amine trically opposite 
the molecular beam. For analyzing the direct beam, the 
pulse travels directly to a flow-through type ionization 
detector as shown by Position A of Fig. 3. 1.2. For stay- 
time measurements, the pulse first travels to a target and 
then to the detector as shown by Position B. The signal 
from the ionization detector, proportional to instantaneous 
number density within its active zone, is amplified and then 
fed into a signal averager. The averager divides the signal 
into 100 consecutive channels of equal time width and 
averages the signal over many pulses to improve the signal- 
to-noise ratio. The averaged signal is then recorded in 
analog and/or digital form. 

The time evolution of the detector signal depends on 
the initial shape of the pulse, the speed distribution of 
beam molecules, the experimental geometry, and for stay-time 
measurements, the distribution of molecular stay times on 
the surface and the speed distribution of molecules leaving 
the surface. The dependence of the signal shape on these 
factors is analyzed in Chapter 4, 

3.3 Targets and Temperature Control 

Three targets were employed. One was high purity 
(99.999$) copper and two were high purity (99.995$) nickel. 
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Each was 0.64 cm diam x 0.16 cm thick and was mounted on a 
copper holder as shown in Fig. 3.3.1. The copper holder was 
attached to a liquid-nitrogen-filled reservoir by means of a 
boron nitride spacer. This spacer electrically insulated 
the target (to allow measurement of the ion current during 
ion bombardment) and still provided adequate thermal conduc- 
tance for cooling the target. A tungsten filament was 
inserted inside the holder to allow heating of the target 
and was electrically insulated from the holder by a ceramic 
tube. A chromel-constantan thermocouple was mounted on the 
holder just below the target for measuring the target 
temperature. The output of this thermocouple was recorded 
on a continuous recording potentiometer and was also used in 
a feedback loop to control the voltage across the heater 
filament and thereby control the target temperature. With 
this arrangement, the target temperature could be varied 
from 90°K to 450°K and could be maintained constant within 
+0.5°K. 

3.4 Target Preparation 

All targets were first polished to a mirror finish, 
following standard metallurgical practice, finishing with 
No. 600 diamond paste. The Cu target and one Ni target 
(hereafter referred to as Ni A) were then cleaned by rinsing 
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Figure 3.3.1 Target Mount 
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in acetone followed by distilled water. After mounting the 
target in the vacuum chamber and evacuation of the chamber, 
the target was heated to about k 50 °K and held there for 3 
hours. The target was then cooled and stay-time experiments 
were performed. The second Ni target (hereafter referred to 
as Ni B) was heated in a vacuum oven to 1300°K for 30 
minutes, removed, and then subjected to a series of cleaning 
steps known as the "Diversey" process. Although the detailed 
compositions of the cleaning agents used in this process are 
proprietary to the Diversey Company, they consist basically 
of (a) a degreaser, (b) an oxidizer, (c) an oxide remover, 

(d) an etchant and finally (e) a distilled water rinse. 

This process has been used successfully in other experiments ^ 27 ^ 
on nickel to remove all but a few atomic layers of oxide 
from the surface. The Ni B target was then subjected to 
vacuum heating and stay-time experiments similar to those 
for the Cu and Ni A targets. After a number of experiments 
with Ni A, it was removed and cleaned by the Diversey pro- 
cess. 

For clarity, the steps to which each target was 
subjected are listed in Table 3.^.1. In this table, speci- 
fic target preparation procedures are denoted by a shorthand 
notation. An observable difference between Ni A (or Ni A+D) 
and Ni B + D was the larger grain size for the Ni B + D 
resulting from the vacuum firing. An observable difference 
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between Ni A and Ni A + D resulted from the etching part of 
the Diversey process. 

In addition to the above preparations, some of the tar- 
gets were subjected to ion bombardment. The ion gun is 
shown in Fig. 3.4.1. In this gun, high purity (99.995#), 
inert gas (He, Ar, or Xe) is injected into the gun’s ioniza- 
tion region. Positive ions are then accelerated out of this 
region, pass through a series of focusing rings, and emerge 
with an energy of 300 eV. The ion beam strikes the target 
at an angle of incidence of about 45°. From the results of 
preliminary measurements with this gun in a vacuum bell $ar, 
the ion beam size at the target is estimated to be no 
greater than 1.0 cm diameter. Ion currents to the target 
were measured by appropriately connecting the thermocouple 
leads through an electrometer. Using argon gas in the gun, 
ion currents were typically 2 x 10“ 9 amps. The gas injec- 
tion rate was such as to cause the test chamber pressure to 
rise by 2 x 10”^ torr. 
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Figure 3.^.1 



INSULATED BARREL 


- ION BEAM 


TZZZZZZZSZZZZZZZi 


FOCUS 

VOLTAGE 



Ion Gun 


CHAPTER 4 

RELATION OF DETECTOR SIGNAL TO MEAN 
STAY TIME FOR A FIRST-ORDER DESORPTION PROCESS 

The shape of the detector signal for a distribution of 
stay times corresponding to a first-order desorption process 
is examined in this chapter. The effects of finite shutter 
open time and finite molecular transit times from chopper- 
to-target and target-to-detector are included in this 
examination for the geometry and experimental conditions of 
the present study. For the assumed first-order desorption 
process, a method is presented for deducing mean stay times 
from certain parameters of the detector signal. The analy- 
sis is a generalization of that presented by Bailey^ 1 ^^ to 
include the effects of finite shutter open time and incident 
beam speed distribution. 

4.1 Detector Signal for a Triangular Shutter Function and 

First-Order Desorption 

Consider the case of unsteady adsorption followed by 
desorption at a rate proportional to the surface coverage 
(first-order desorption). The total rate of change in 
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surface coverage is then 



dt 


p a p i 



T 


(4.1.1) 


If the adsorption rate, given by P a F^(t), is known, then the 
solution of Eq. (4.1.1) is 


<r a (t) = exp(-t/r) / P a F^(t* Jexptf/rJdt* 

— oo (4.1.2) 


The stay time r is the mean stay time for the ensemble 
P a F^(t') dt' . The desorption flux F d (t)» of course, is 
simply 


F d (t) = * a (t) /t 


(4.1.3) 


The effects of shutter open time and transit times can 
now be added to determine the expected shape of the detector 
signal. The shutter time and transit time from chopper-to- 
target are contained implicitly in the incident pulse F^t). 
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The shape of this pulse is dependent upon the speed distri- 
bution for the incident beam and upon the shutter function 
g(t), (The shutter function describes the shape of the 
pulse as it emerges from the chopper slot and depends on 
the rotational speed of the chopper and the beam-chopper- 
target geometry, ) In the present experiments, the effective 
open time of the shutter is much larger than the mean tran- 
sit time from the chopper to target. In Appendix III, it is 
shown that in this limit 

F i (t) ~ g(t-t Ql ) (4.1.4) 

where * OJ is the most probable time of flight from the 
chopper to the target. If g(t) is approximated by a symme- 
tric triangle of basewidth 2t g , then 

l+(t-t 0l )/t s for -t s < t-t 0l < 0 
l-(t-t 0l )/t g for 0 < t-t 0l < t a 
0 otherwise 

(*. 1 . 5 ) 

where F Q is an appropriate proportionality factor (approxi- 
mately equal to the steady-state beam flux striking the 
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target if the chopper were removed). 

If the desorbing molecules have a Maxwellian speed 
distribution corresponding to the surface temperature, it 
can be shown^ 1 ^ that the detector signal S(t) will have 
the form 


S(t) = S . 


/ 


a aS v > 


-*.**01 


fey exp fe(s 


dt* 
(4.1.6) 


where 0 g =* sfcT/ 2kT s , x 2 is the target-detector distance and 
S o is 811 a PP ro Priate proportionality factor (since only the 
time-wise shape of the signal is of interest, the propor- 
tionality factor is unimportant). The time p^x^/jz is the 
most probable transit time from target to detector. In the 
present experiments, t g and P q x 2 were about the same magni- 
tude and both were greater than t Ql by a factor of 15 to 20. 

4.2 Relationship Between t and Parameters of the Predicted 

Signal Shape 

The shape of the detector signal is now obtained by 
combining Eqs. (4.1.2), (4.1,5) and (4.1.6). At this point, 
it will be further assumed that P & * 1. Substituting F i (t) 
from Eq, (4,1.5) into Eq. (4,1,2) and integrating, one 
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obtains a (t). (The exact result of this, integration is 
given in Appendix III along with some typical shapes for 
<7 (t) calculated from this expression.) Then substituting 
this result for a A t) into Eq. (4.1.6) yields an integral 

a 

expression for S(t). 

Some typical predicted signals, obtained by numerical 
integration of the above described expression of S(t), are 
shown in Pig. 4,2.1. All times shown in this figure are 
non-dimensionalized by the quantity 6 ®v^ o1 . The values of 
t /b and /?x«/d are typical for the experiments. The time 
t at which S(t ) is a maximum and the time t 0 at which 
S(t) = Js(t m ) are used to characterize the signal shape. 

For non-zero stay time, a "peak shift" 


^ t ~ t* — t _ 
s m m,T-o 

and a "time constant" 


(4.2.1) 


dt 


T ^ = . 


d[ln S(t)] 


t=t. 


(4.2.2) 


are more direct measures of the stay time. For small 

T ( < t m T3 . 0 « 0 g x 2 ), ** shown * as might have been 

expected, that the detector signal is simply delayed by t, 

and therefore the peak shift At a is equal to the mean stay 

s 
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time. For T»t and for t »t the detector 

nifV -o m f T=o 

signal simply decays as exp(-t/r) and the time constant r z 
equals the mean stay time. The behavior of dt and tu as 

S 2 

functions of t is illustrated in Fig. 4.2.2 for typical 
values of t g and /? g x 2 (all times are again normalized by $). 
This figure shows that, for intermediate values of t, 
neither At g nor is a good measure of the stay time. By 
observing the slope of 4t g versus T, it is also seen that 
even though the shutter open time is finite or even large 
compared to r(i.e., t » t) » the peak shift At is still 

® S 

equal to the stay time for small t . 

In general, the predicted expression for S(t) and the 

quantities At g and are functions of the parameters t , 

s 

* 01 ' ^s x 2 and T ' For a Siven detector signal, all except t 
will be known from other experimental information. Thus, if 
the derived model for S(t) were sufficiently accurate, t 
could be determined either by some sort of curve fit of the 
predicted function to the measured signal shape or by com- 
parisons of predicted to measured At g or t 2 . Of course, the 
assumptions made in deriving this model are only approxima- 
tions, and some disagreement between measured and predicted 
signal behavior is expected. At least in principle, the 
validity of these assumptions can be assessed by comparisons 
of the predicted behavior to that observed experimentally. 








CHAPTER 5 

EXPERIMENTAL PROCEDURES 


The discussion of the experimental program will be 
divided into two parts i (a) calibrations and (b) stay-time 
experiments. Detector calibration and characterization of 
the molecular beam and TOP detection system are described in 
Appendix II. The stay-time measurements and the variables 
involved are described in Sec. 5.1. The techniques for 
extracting stay times and related parameters from the 
recorded detector signals are described in Sec. 5.2. An 
analysis of the uncertainties associated with the measure- 
ments and the data reduction techniques is given in Sec. 5.3. 

5.1 Stay-Time Experiments 

The primary stay-time experiments are listed in Table 

5.1.1 along with the experimental conditions for each run. 

The target designations are those given in Table 3.4.1. The 
ranges of target temperatures listed are those for which 
significant peak shifts were observed. For calibration pur- 
poses (see Sec. 5.2), detector signals were also recorded at 
a number of temperatures between the upper limits listed and 

room temperature. Each recorded signal represents an aver— 

4 

age over about 10 repetitive pulses (a period of about 2 
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TABLE 5.1.1 


PARAMETERS FOR THE STAY-TIME EXPERIMENTS 


RUN 

DESIGNATION 

GAS 

SPECIES 

TARGET 

— 

TARGET 

TEMPERATURE 

°K 

ION 

BOMBARDMENT 

SPECIES 

ION 

DOSAGE 
PRIOR TO 
RUN, C 

ION 

CURRENT 
DURING 
RUN, A 

Xe — —Cu 

Xe 

Cu 

117—102 

None 

0 

0 

Xe -Ni A 

Xe 

Ni A 

123—105 

None 

0 

0 

Xe -Ni A + D 

Xe 

Ni A + D 

117—105 

None 

0 

0 

Xe »N1 A + D + He + 

Xe 

Ni A + D 

111— 99 

He 

5 X 10' 6 

2 X 10“ 9 

Xe ► Ni A + D + Ar + 

Xe 

Ni A + D 

117 — ^ 99 

Ar 

5 x 10' 6 

3 x lo" 9 

Xe - Ni A + D + Xe + 

Xe 

Ni A + D 

105—99 

Xe 

, ‘ 6 
5 x 10 

-9 

2 x 10 

Xe — Ni B + D 

Xe 

Ni B + D 

123—102 

None 

0 

0 

Xe — Ni B + D + At 4 

Xe 

Ni B + D 

114—95 

Ar 

5 x l(f 6 

-9 

2 x 10 

Kr ^Ni A + D 

Kr 

! Ni A + D 

99 — 92 

None 

0 

0 

COg— Ni A + D 

co 2 

Ni A + D 

123—108 

None 

0 

0 
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minutes) during which time the target temperature was held 
constant. 

In addition to the experiments listed in Table 5.1.1* 
detector signals were recorded for Xe-*NiA at T g = 111°K at 
various source pressures between 1 to 10 torr. Detector 
signals were also recorded as a function of exposure of the 
surface to various gases for Xe -► NiA+D+Ar + at T g * 108°K. 
Other experimental conditions were the same for all experi- 
ments listed in Table 5.1.1 and are listed in Table 5.1.2, 

5.2 Data Reduction 

Peak Shifts . In order to Hind the previously defined 
peak shift for a measured detector signal, it is necessary 
to determine the peak location t for that signal and the 

Ul 

time t corresponding to the temperature at which the 

in i / ™o 

signal was measured. The time t m for each recorded signal 
was determined from a least-squares-polynomial fit to the 
recorded signal in the neighborhood of the apparent maximum. 
To establish the quantity 0 , the expected variation in 

this quantity with surface temperature was used. Experimen- 
tally, it has been found that, for T at or below room 

s 

temperature, most gases thermally accommodate themselves to 
the solid before leaving the surface. Thus, in the range of 
T g where t is negligible and with negligible shutter open 



TABLE 5.1.2 


FIXED CONDITIONS FOR THE EXPERIMENTS 


SOURCE PRESSURE, p 

’ o 

SOURCE TEMPERATURE, T 

’ o 

SHUTTER OPEN TIME, t 

’ s 

SHUTTER -TO-TARGET DISTANCE, Xj 

TARGET-TO-DETECTOR DISTANCE, x g 

ANGLE BETWEEN INCIDENT BEAM AND 
TARGET NORMAL 


10 t 0. 1 torr 
295 ± 3 °K 
2.4 ± 0.4 X 10' 
0. 32 ± 0.05 cm 
4.36 + 0.1 cm 
45 t 5 0 


ANGLE OF REFLECTION INTO DETECTOR 45 t 5 ° 

MOST PROBABLE TRANSIT TIME FROM 1.1 + 0 2 X 10~ 5 

SHUTTER TO TARGET, 0 x ly </2 
(Calculated for Xe) 0 


sec 


sec 


-p- 

tv> 



43 


time, one would expect 

Vo * A + sA/r B (5.2.1) 

where A and B are constants. However, for finite shutter 
open time, it is not obvious that t predicted from Eq. 

m f T=“0 ^ 

(4.1.6) will show precisely the above behavior. Therefore, 
t n was measured for three gases, Xe, Kr and C0 2 , over a 
range of T g where r was expected to be negligible. Corre- 
sponding values of t were also determined for detector 

uj f * “ O 

signals predicted by Eq. (4,1.6). Both sets of results are 
shown in Pig, 5.2.1 and t TaQ is seen to vary linearly with 
1/VV The indicated quantitative agreement between 
measured and predicted values is an artifact to some extent. 
The measured values had a fixed uncertainty due to a time 
displacement between g(t) and the chopper photopulse. (This 
uncertainty was virtually unimportant in determining peak 
shifts, since the displacement was presumedly constant for a 
given data run.) Thus, the measured t^* s shown in Pig. 

5.2.1 have been adjusted such that the respective best-fit 
straight lines through measured and predicted values have 
the same intercept A at 1/i/r^ = 0, Nevertheless, the temper- 
ature dependence of Eq. (5.2.1) is quite close to that 
observed and therefore provided satisfactorily accurate 
estimates of t m T=0 . Values of A and B were determined from 




Figure 5.2.1 Compa: 
Most ] 
Targe - 
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data for each run similar to that shown in Pig. 5.2.1. 

Time Constants . The time constant, t 2 , was determined 
by fitting the expression 

S(t) = Ce~ t/V 2 (5.2.2) 

to the detector signal in the neighborhood of t = t 2 « The 
above form, equivalent to the definition given in Eq. 
(4.2.2), was more convenient to use in computerized data 
processing. 

Mean Stay Time . The stay-time t was estimated from 
dt using the detector signal model given in Chapter 4. 
Given values of t g , t Ql , and 0 g x 2 , an initial value for r 
was assumed and used to calculate S(t) from Eq. (4.1.6). 

The value of dt g for this S(t) was determined and compared 
to the experimental value. Iteration on t was continued 
until the calculated and experimental values of dt g agreed. 
The quantities t Ql and /9 g x 2 used in the above were calcu- 
lated using the measured source and target temperatures, 
respectively, and the flight lengths x x and x 2 . The effec- 
tive shutter open time t was determined from a separate 
series of measurements in which detector signals were 
recorded as a function of chopper speed. In the limit of 
large shutter open time (t g >>£ g x 2 and r), it can be 
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( 26 ) 

shown' ' that the detector signal shape will approach that 
of the shutter function g(t). For a triangular shutter 
function, t g is just the width at half-maximum and should be 
linearly proportional to the chopper period. The observed 
variation of the detector signal width at #S(t m ) was also 
found to be linear with chopper period at slow chopper 
speeds. Thus, this variation was used to estimate the 
effective open time at higher chopper speeds of the trian- 
gular shutter function used in the stay -time estimations. 

Binding Energy and Pre exponential Factor . If the stay 

time obeyed a Frenkel-type relation and if the first-order 

detector signal model were perfectly accurate, then the 

measured dependence of the above determined t with T could 

8 

be used to obtain the binding energy E Q and pre exponential 

f 

factor t 0 . However, as will be seen in Chapter 6, disagree- 
ment between measured and predicted detector signals gave 
reason to question some of the model assumptions. As a 
result of this disagreement, the stay times calculated from 
the peak shifts are considerably in error at the lower sur- 
face temperatures. Therefore, the qualitative behavior of 
the factors E Q and r Q was determined directly from the peak 
shift data assuming 


' lt s a T o •*P(V k V 


( 5 . 2 . 3 ) 



The parameters, E Q and T 0 , were estimated for each data set 
by the following procedures! 

(1) Direct curve-fit of Eq. (5.2.3) to peak shift 
data* 

(2) Calculate E 0 for each data point based on 

mobile-adsorption prediction, Eq. (2.1,11), 

assuming t « At » 
s 

(3) Calculate E Q for each data point based on 

local! zed-adsorption prediction, Eq. (2.1.12), 

again assuming t » dt . 

s 

The quantities <? D and <7^ required in procedures (2) and (3) 
were taken from Refs. (28) and (29) respectively. 

5.3 Uncertainty Analysis 

Uncertainties in the results are concerned mainly with 
the following quantities. First are the uncertainties in 
the experimental parameters, particularly the surface temper- 
ature. Second are the uncertainties in stay- time parameters 
deduced directly from the recorded detector signals* these 
errors are introduced both by the actual measurement process 
and the associated data reduction techniques. Next are the 
undertaint ies in stay times, binding energies and preexponen- 
tial factors that are derived from the measured stay-time 
parameters and experimental parameters. These latter errors 
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depend on the parameter uncertainties but are expected to be 
mainly due to approximations in the detector-signal model 
and adsorption model employed in deducing said quantities. 
Finally are uncertainties in the nature of the surface which 
make data interpretation more difficult. 

Uncertainties in experimental conditions . The primary 
experimental variable which affects the stay time is the 
target temperature. Systematic error in reported values of 
the absolute target temperature is less than 3°K (based on 
calibration of the target thermocouple at the boiling point 
of liquid nitrogen). However, from the reproducibility of 
the indicated temperature at different constant heat inputs, 
errors in temperature differences are estimated to be less 
than +0. 5°K. Target temperature fluctuations during the 
averaging of a given detector signal were typically +0.5°K. 
Other parameter uncertainties which affect the detector 
signal are listed in Table 5.1.2. 

Uncertainties in measured stay-time parameters . Errors 
in dt g and t 2 result due to timing errors in the signal 
measurement and due to extraction of these parameters from 
the recorded signals. The time response of the detector and 
its signal processing instrumentation was determined from 
time-of-flight measurements with a siggle-orifice molecular 
beam source (see Appendix IX). Based on these measurements. 
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timing errors due to signal shape distortion by the detection 
system are estimated to be less than 10 J sec. The accuracy 
of the data reduction for dt s and Tg depends primarily on 
the signal-to-noise (S/N) ratio of the recorded signal. 

Since the S/N generally decreased with decreasing tempera- 
ture, the absolute uncertainty in these parameters is larger 
at the lower temperatures. The uncertainty in determining 
t for a given recorded signal varied typically from about 

±10“^ to ±10“ 4 sec (or +3 to ±6# of measured t ). Based 

in 

strictly on the observed linearity (and deviations from 

linearity) of t with 1//T7 at temperatures where t was neg- 
in s 

ligible, the uncertainty in predicting t m T=0 at lower 
temperatures by Eq. (5.2,1) was about ± 10 “^ S ec (or ±2% of 

measured t ) . These uncertainties in t and t „ rep- 

resent the limits of 95# confidence in the data reduction 
techniques employed. The resulting uncertainties in the data 
reduction for the peak shift, estimated from the rms sum of 
the above, varied from about ±2 x 10”^ to ±10“ 4 sec. The 95# 
confidence limits in determining were more sensitive to 
the S/N ratio and varied from about ±4 x 10“^ sec at the 
higher temperatures to about ±10'^ sec at the lowest temper- 
atures for these experiments (corresponding to ±10 to ±25# 
of measured t- 2 ). 

f 

Uncertainties in t, E q and T' , The accuracy of the 
stay time estimated from At e is primarily dependent on the 

o 
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validity of the detector-signal model. The approximations 

* 

employed for the incident and desorbed beam speed distribu- 
tions are expected to be sufficiently accurate for determining 
r in this way. The accuracy of the triangular shutter func- 
tion approximation and the validity of the first-order 
desorption model will be discussed in the next chapter. 
Geometrical effects not considered in the model, such as the 
finite size of the target and detector aperture, are esti- 
mated to have negligible effect on the signal shape and thus 
negligible effect on t. 

The uncertainties in the parameters, t ' and E . are 

o o' 

dependent on the validity of the adsorption model and will 
be discussed in Chapter 6 along with those results. 

Uncertainties in surface conditions . The composition 
and physical structure of the surfaces of the various tar- 
gets is subject to considerable uncertainty. The coverage 
due to gases (other than the beam species) adsorbing from 
the test chamber background was probably significant parti- 
cularly at the lower temperatures. Nevertheless, the 
reproducibility of these conditions is verified to some 
extent at least by the reproducibility of the experimental 
results presented in Chapter 6. 



CHAPTER 6 
RESULTS 


6.1 Comparison Between Measured and Predicted Detector 

Signal Shapes 

Three typical recorded detector signals are shown in 
Fig, 6.1.1, For part (a) of this figure, the detector sig- 
nal was recorded with the target at room temperature where 
the stay time should be very small. The signal shape pre- 
dicted by Eq. (4,1.6) for T = 0 is shown here for comparison. 
The agreement between measured and predicted signals is good 
up to t«t n , but for t >t n) , the measured signal decays at a 
rate smaller than predicted. The disagreement at large 
times may be due in large part to errors in the approximate 
shutter function used in generating the theoretical curve. 

The assumed shutter function is symmetric about its midpoint. 
Calculations by 01ander (30) indicate that the symmetry of 
the angular flux distribution and therefore the symmetry of 
the shutter function are sensitive to source-collimator 
alignment. Some additional calculations of signal shapes 
were performed using an asymmetrical triangular shutter 
function. By varying the degree of asymmetry in this func- 
tion, a shutter function was found which predicted the 
measured signal shape up to twtg. Comparing this 
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Figure 6,1.1 Comparison of Measured and Predicted 
Detector Signals, Xe-*-Ni A 
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asymmetry to predictions of Olander, it is estimated that 
either a source displacement of about 0.5 source -diameters 
from the nominal beam axis or a source -exit -plane tilt of 6° 
could explain the discrepancies in Pig. 6.1.1(a). Such 
errors are within the misalignment uncertainties for the 
present beam geometry. 

In principle, any measured signal for t = 0 could be 

used to solve for the exact shutter function. However, for 

small t , it is the shift in the peak with increasing stay 

time that is of interest. The observed variation in t with 

m 

target temperature for negligible t is that expected. Thus, 
to first approximation, the disagreement between theory and 
experiment at large t should introduce negligible error in 
the determination of t from peak-shift data. Of course, 
estimation of t from the observed decay rate of the signal 
at large t (i.e., t = tg) will be affected. Even so, as t 
becomes large compared to the discrepancy in time, the error 
introduced in t 2 will approach zero. 

A second typical signal, measured at a target tempera- 
ture where the peak shift is finite due to stay time, is 
shown in Pig. 6.1.1(b). Also shown here are the predicted 
signal for t * 0 and the predicted signal for that value of 
r which fits the measured signal at t . To first order, the 
predicted (t ^ 0) curve is delayed approximately by t up to 
t «t ffl , but beyond t ffl the time delay becomes increasingly 
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greater. The predicted {t/ 0) curve agrees with the measured 

curve up to t«t . However, for t >t„, the measured curve 

in in 

decays at a lower rate than the predicted (t ^ 0) curve. 

This discrepancy at large t is similar to that observed at 
room temperature for negligible stay time. 

At even lower target temperatures (larger stay times) 
the measured detector signals take on a shape that was alto- 
gether unexpected. An example of this shape is shown in 
Fig. 6.1.1(c). It is seen that the signal exhibits a "nega- 
tive" dip at small times. To understand a negative signal, 
it must be recalled that the detector's ion-collector cur- 
rent includes a component due to background gas in the 
chamber as well as a component due to the pulsed molecular 
beam. During detector signal processing, only the rela- 
tively fast, time-varying components are amplified. That is, 
the low-frequency (below about 10 Hi) components are inten- 
tionally filtered out. Moreover, signal averaging eliminates 
any background fluctuation which does not occur repetitively 
at the chopper frequency. ThuB, under normal conditions, 
the background component does not appear in the recorded 
signal. However, the background gas is composed, in large 
part, of the same gas as in the beam and is steadily adsorb- 
ing and desorbing along with the beam pulses, Sinoe it is 
impossible for the beam pulses to give rise to the negative 
dips, it is most likely that the beam pulse in some way 



temporarily reduces the rate of desorption of background 
molecules from the target. 

Signals exhibiting the negative dip obviously cannot be 
described by the detector- signal model for first-order desorp- 
tion. Still, it was observed that both dt g and -r 2 for this 
type of signal increased with decreasing surface temperature 
in a reproducible manner for a given data run. In addition, 
the signal at long times (t~t 2 ) seems to decay exponentially 
as for a first-order process at least to the extent that can 
be determined by curve -fit techniques. Thus, dt_ and t 0 
appear to be indicative of at least the qualitative behavior 
of the stay time even for this type of signal. The behavior 
of signals exhibiting the negative dip is discussed further 
in Sec. 6.5. 

6.2 Stay-Time Parameters 

Stay-time parameters dt g and v 2 deduced from the set of 
detector signals from which the ones in Pig, 6.1.1 were 
taken are shown in Pig. 6.2.1. Also shown are values of t 
derived from the dt g data using procedures described in 
Chapter 5. The error flags shown represent the limits of 
95# confidence for determining dt and 'r. from the recorded 

o c. 

signals by the data reduction techniques described in Chap- 
ter 5. (Note that the scatter in the plotted data generally 
falls well within these limits.) If the stay time varied as 




Figure 6.2.1 Typical Variation of Measured Stay-Time 
Parameters with Target Temperature, 
Xe-*Ni A, Error Flags Denote 95 % 
Confidence Intervals in Data Reduction 
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exp(E Q /kT g ) and the negative dip did not occur, dt g and x 2 

should exhibit variations with T” 1 similar to the variation 

s 

with r shown in Pig. 4,2.2. For dt < lO”^ sec, this is the 

s 

-4 

case. However, for dt > 10 sec, x derived from dt shows 

s s 

a rapid increase as T* 1 increases and eventually exceeds X 0 . 
Since in theory one expects x 2 /x> 1, these x values are 
clearly too high at large T' 1 . 

3 

The linearity of lndt_ with T" 1 noted here was also 

s s 

observed in all other experimental data. At small x, the 
slope of lndt g must be the slope of lnr with T" 1 . At large 
r, the x 2 data indicate that the actual stay time is no more 
than a factor of 4 greater than dt 0 . Thus, it seems that 
the dt g data are better representations of the variation of 
the actual stay time over the temperature range than one 
should expect. For this reason, only the peak-shift (dt ) 

9 

data will be considered in the following presentation of the 
effects of other parameter variations. However, x 2 and x 
data comparable to those shown here have been determined for 
each recorded signal and are tabulated in Appendix IV. 
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6.3 Results of Parameter Variations 

6.3.1 Effect of Target Material for Xe Adsorption 

Peak shift data for Xe adsorption on Cu and NiA are 
plotted in Pig. 6.3.1. It is seen that the peak shifts 
obtained with NiA are about twice those for Cu at any given 
T. Since both targets are probably covered to some extent 
by other physisorbed or chemisorbed species, the differences 
may be due as much to differences in adsorbates as to differ- 
ences in the target bulk material. At least one other run 
for each target reproduced the data shown within approxi- 
mately the symbol size shown in this figure. Thus, the 
observed peak shift differences are outside the repeatability 
of the data. 

6.3.2 Effect of Pre-Vacuum Surface Preparation for Xe 
Adsorption 

The effects of the previously described pre-vacuum sur- 
face preparations for the nickel samples are illustrated in 
Pig, 6.3.2. The variation of dt with T" 1 for Xe — NiA and 

8 S 

Xe-*NiB+D are almost identical. The best-fit line for 
Xe -* NiA+D data intersects the best-fit lines for the other 
data but has a larger slope. This may be related to the 
known roughening of the NiA surface by the etchant in the 









1/T S , °K’ 1 


Figure 6.3.2 Variation of Peak Shift with Target 
Temperature for Different Pre -Vacuum 
Surface Preparations 
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Diversey process. (While the NiB+D was also etched, its 
grain size prior to etching was much larger than NiA due to 
grain growth during vacuum firing. Therefore, both NiA and 
NiB+D were smoother to the eye than NiA+D.) The higher 
slope for NiA+D implies a higher heat of adsorption and this 
is in qualitative agreement with other work (see for example 
the discussion of surface non-uniformity by Young and 
Crowell (1) ). 

6.3.3 Effect of Argon Ion Bombardment During Xe Adsorption 

The effects on the peak shifts of argon ion bombardment 

of the target are shown in Fig. 6 . 3 . 3 (a) for Xe^NiA+D and 

in Fig. 6. 3. 3(b) for Xe-*NiB+D, At fixed T , ion bombardment 

s 

leads to a factor of 10 to 20 decrease in the peak shift. 

This result is somewhat surprising, since ion bombardment is 
usually thought of as a means for creating a cleaner surface 
and therefore presumably one with higher binding energy for 
adsorption. 

The surface changes that are most likely to occur dur- 
ing ion bombardment are due tot 

(1) sputtering of adsorbate 1 

(2) sputtering of bulk material* 

(3) implantation of neutralized ion species. 

For the present experiments, the ion flux was much less than the 
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flux due to either the bean or the background gas. Thus, 

ion bombardment should not significantly alter adsorbate 

coverage from these sources. Although the ion flux was low, 

the surface was exposed to this beam for a significant time. 

Based on a sputtering probability of unity^ 1 ^ for 300 eV 

Ar+, it estimated that about 0.1 of the surface atoms 

should have been sputtered. However, such sputtering should, 

if anything, increase the effective binding energy and 

therefore increase the stay time. Thus, the most likely 

explanation for the observed decrease in peak shift is that 

the argon ions bury themselves near the surface and reduce 

the binding energy between the surface and the Xe atoms. 

For 300 eV argon ions, Johnson w/ estimates the average 

penetration distance in nickel to be about 4X. The high 

( 33 ) 

desorption energies ' expected for atoms (neutralised 
ions) buried even one atomic layer from the surface would 
result in a negligible desorption probability at these tem- 
peratures, The degree to which these buried atoms could 
affect the binding energy depends on their number density. 

For a total ion dose of 5 x 10~^C P the argon surface cover- 
age could be at most about lO 1 ^ atoms/cm 2 , about 0.1 
monolayer. If the net adsorbate -solid attraction were 
simply a sum over two-particle attractions and assuming that 
the attractive force between Xe - Ar is much less than 
between Xe - Ni, a coverage of 0,1 monolayers of Ar might at 
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most give rise to a 10# decrease in binding energy for Xe. 
Because of the exponential dependence of t on E Q , a 10# 
decrease in E Q can easily result in a factor of 10 decrease 
in t. Thus, the observed effect of ion bombardment could 
easily be attributed to ion burial near the surface, 

6.3.4 Effect of Ion Bombarding Species for Xe Adsorption 

Following Ar bombardment, the NiA+D target was bom- 
barded with He and then Xe ions. In Fig, 6.3.4, peak shift 
data measured during ion bombardment with each species are 
compared to those observed prior to ion bombardment. He and 
Ar ions are seen to produce comparable reductions in the 
peak shift, while Xe ions produce an even further decrease. 
It was observed that, after the ion beam was shut off 

and with T_ fixed, the peak shift would increase with time, 
s 

The rate at which this increase occurred was quite low» even 
after several days, the peak shift would not return com* 
pletely to its value prior to ion bombardment. This 
indicates that a relatively permanent change in the surface 
occurred during ion bombardment. Thus, the chronological 
order of the tests could be very important. Time dependence 
of the measured stay time under various experimental condi- 
tions is discussed further in Sec. 6.3.7. 









69 


6.3.5 Effect of Incident Flux Variation for Xe Adsorption 

The effect on the peak shift of varying the incident 
flux, accomplished by varying the source pressure (flux 
proportional to N /p^ 2 '^)« is shown in Fig. 6.3.5. The peak 
shift is seen to increase by about 5 0% between p Q - 1 torr 
and p^ = 3 torr and to reach a broad maximum near p^ = 6 torr. 
In the absence of adsorbate -adsorbate interactions, 7- should 
be independent of p c . However, the appearance of the pre- 
viously noted negative dip in the detector signal indicates 
that such interactions are occurring. An estimate of the 
total surface coverage of Xe from both the beam and back- 
ground was made using the incident flux of J.k x 10*-* 

-2 -i 

atoms • cm .sec given in Chapter 3 for p Q = 10 torr and 

the stay time estimated from the peak-shift data in 

Fig. 6.3.5. The value obtained for this coverage was about 
11 -2 -3 

5 x 10 atoms .cm or about 10 J monolayers of Xe. Such 
low coverage would not be expected to give significant 
Xe - Xe interaction. Of course, ads or bate -ads or bate inter- 
actions may have occurred between Xe and some other species 
adsorbed from the background, but again due to the low 
coverage of Xe, the peak shift should be independent of p o * 
Thus, in the absence of additional information about the 
state of the surface, the observed dependence of 7 - on p 

o 

cannot be interpreted. 
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6.3.6 Stay Times for C0 2 and Kr Adsorption 

Peak shift data for C0 2 and Kr adsorbing on NiA are 
compared to those obtained for Xe in Fig. 6.3.6. The peak 
shifts for C0 2 are comparable in magnitude to those for Xe 
but exhibit a smaller slope versus T" 1 . The temperatures at 
which Kr shows comparable stay times are about 20# lower 
than those for Xe, This would be expected on the basis of 
known differences in atom-atom attractive well depths for Xe 

( 34 ) 

and Kr w , While the molecule -molecule well depth for 
(34) 

C0 2 w falls between those for Xe and Kr, it is not sur- 
prising that stay times for a polyatomic gas do not scale 
exactly with this well-depth. 

6.3.7 Variation of Peak Shift with Time 

Bailey * 15 * observed that after cooling his surface to 
a given temperature, the stay time decreased with elapsed 
time. This decrease was attributed to the adsorption of 
gases present in his chamber. Since the present vacuum 
system and pressures are similar to those in Bailey’s work, 
the time dependence of the peak shift was investigated in 
this work also. 

First, the peak shift was measured as a function of 
time both with and without simultaneous ion bombardment. 

These data are shown in Fig. 6.3.7. Also shown in this 




1/T S , °K' 1 


Figure 6.3.6 Variation of Peak Shift with Target 
Temperature for Different Incident 
Gas Species, Ni A+D Target 




Figure 6.3.7 Variation of Peak Shift with Time 
During Argon Ion Bombardment for 
Different Initial Ion Dosages, 
Xe-Ni A+D, T s = 108°K 
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figure are data obtained for different accumulated ion doses 
at time = 0. (The rate prior to t = 0 for these data was 
that shown for t>0. This accumulated dose is the total 
charge due to ions received by the sample immediately prior 
to recording the dataj it is not the total dose the target 
has received since initial installation. ) 

Por no ion bombardment during the recording of the peak 
shift data (circles), the peak shift shows a gradual 
increase with time. The target received an accumulated dose 
of about 5 x 10 on the day previous to this run, and this 

increase with time is similar to that noted in Sec, 6.3.4 
upon cessation of ion bombardment. However, similar 
increases in dt 0 with time were observed for targets having 
received no previous ion bombardment. 

For ion bombardment beginning at t = 0 (squares), the 
peak shift first decreases sharply with time, levels off, 
reaches a second maximum and decreases toward some asymptotic 
limit. The sharp initial decrease is probably due to a 
larger ion gun current known to occur during the first few 
minutes of operation. Por an initial ion dose of 2 x 10“^c 
(diamonds), dt g shows a gradual decrease with time, leveling 
off after about 26 minutes. After an initial dose of about 
5 x 10 C at t « 0 (triangles), dt g remains essentially 
constant. Assuming the previously discussed ion -burial 
theory is correct, the data in this figure suggest that the 
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surface becomes saturated after an ion dose of about 

-9 

5 x 10 C and that bombardment at a rate of 2 x 10 A can 
maintain saturation. 

Time variations of the peak shift following cessation 
of ion bombardment are shown in Pig. 6,3.8. For each set of 
data, the ion dose prior to time zero was at least 5 x 10“^ C 
and ion bombardment was stopped at time zero. For Fig. 
6.3.8(a), the sample was exposed only to the residual back- 
ground gas. The background was at a pressure of 3 to k x 
10” 7 torr and consisted mostly of Xe from the molecular 
beam. Other gases, primarily N 2 , 0 2# H^O and C0 2 , as deter- 
mined by a mass spectrometer, contributed a partial 
pressure of about 2 x 10” torr. The effect of this expo- 
sure is seen to be a gradual increase in the peak shift with 
time which continued to very long times. For the data in 
Fig. 6. 3. 8(b) -(d), the partial pressures of C0 2 , H 2 0, and Xe 
in turn were intentionally increased during the measurement 
sequence. For comparison, the data from Fig. 6.3.8(a) are 
shown as a solid line in these latter figures. If the 
increase in peak shift with time was due to adsorption of 
these gases, an increase in the partial pressure of these 
gases should alter the rate at which At increased. Both 
C0 2 and HgO show a possible effect. After about 20 minutes 
exposure to C0 2 , the rate at which At is increasing 
seems to increase (based only on the last data point). 




Time , min 

(a) Background only, p Xe =3x 10“ 7 torr, P other = 2 x 10“ 8 torr 

Figure 6.3.8 Variation of Peak Shift with Time after Ion Bombardment 
for Different Background Pressures and Species, 

Xe -*Ni A+D, T = 108°K 
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After about 20 minutes exposure to HgO, the peak shift 
begins to level off with further exposure time. On the 
other hand, additional Xe exposure, even at much higher 
pressures, shows no significant change in the rate of 4t g - 
increase. Thus, while C0 2 and HgO adsorption may contribute 
to this change in At with time, it seems that some other 
factor is the driving force. It may be adsorption of some 
other species in the chamber or possibly the slow release 
of buried argon from the surface. 

6.4 Binding Energy Results 

» 

The slopes and intercepts, E 0 and v 0 , obtained by a 
least-squares fit of Eq, (5.2.3) to each set of peak-shift 
data are listed in Table 6.4.1. The 95 % confidence inter- 
vals for each curve fit are listed in parentheses in terms of 

decades for t' and in terms of percent of listed value for 
o 

t 

E . The values of t\, are much smaller than typical vibra- 
o o 

-12 _i 4 

tional periods for adsorbed atoms (10 to 10 sec) and 
show large fluctuations (comparable to the confidence 
intervals) from set to set. The values of E Q are physically 
reasonable and fall within the range of values for heats of 
physical adsorption (E^k = 500 to 5000°K). However, the 
values of B Q /k for Xe and Kr on nickel are higher than 
those obtained by Baker and Fox^^ from equilibrium 



TABLE 6.4.1 


PREEXPONENTIAL FACTORS AND BINDING ENERGY TEMPERATURES E Q /k 
OBTAINED BY FITTING PEAK SHIFT DATA TO At = x'expfE /kT ) 

S 0 r 0 s 


DATA SET 

T o » sec 

o 

O 

Ui 

Xe •+ Cu 

6.5xl0~ 18 (10" 19 -I0' 16 ) 

3397 ( + 10%) 

Xe -*• NiA 

6.4xl0~ 17 (10' 17 -10‘ 16 ) 

3213 (± 5%) 

Xe -+ NiA+D 

1.2xlO -20 (10- 24 -10- 16 ) 

4155 (± 24%) 

Xe - NiA+D+He + 

6.3xl0 -17 (10" 18 -10‘ 16 ) 

2916 (± 8%) 

Xe ■+ NiA+D+Ar + 

1.7xl0“ 15 (10' 17 -10‘ 13 ) 

2594 (+ 21%) 

Xe - NiA+D+Xe + 

7.5xI0' 18 (10~ 32 -10 -3 ) 

3010 (+110%) 

Xe - NiB+D 

i.sxio ' 16 oo- 18 -io- 15 > 

3116 (+ 12%) 

Xe - NiB+D+Ar + 

6.3xlO“ 16 (10- 16 -10- 15 ) 

2735 (± 4%) 

Kr -*• NiA+D 

7.0xl0 -20 (10“ 23 -10“ 16 ) 

3278 (± 22%) 

C0 2 -+ NiA+D 

4.2xl0' 18 (10- 19 -10' 16 ) 

3526 (± 8%) 


95% confidence intervals for curve fit are given in parentheses. 
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adsorption measurements who found variations between 2200°K 
and 2?00°K for both Xe and Kr. 

In view of the large uncertainties involved in obtain- 
ing both t q and E Q in the above manner, it was decided to 

r 

estimate t by some other means and then calculate E from 

Q O 

the peak shift data. For this approach, the theoretical 

predictions for mobile and local adsorption were used to 

predict t q and E^k was then computed for each data point. 

The average and standard deviation for the computed E Q /k 

values for each data set are tabulated in Table 6 A. 2 first 

for the mobile adsorption case and then for the localized 

adsorption case. The respective values of r are also 

o 

tabulated there. (Since the r Q for the localized model is 
more sensitive to temperature than t q for the mobile model, 
a range of values for t q is given for the localized case.) 

Theoretical estimates of t' based on the localized- 
adsorption model agree better with the curve-fit estimates 
than do the mobile-model results. It follows that agreement 
For E Q /k is similarly described. Since the observed temper- 
ature dependence for dt g is not grossly different than that 
for the actual mean stay time, it is concluded that the 
localized adsorption model is the better one for these 
experiments. Moreover, the estimates of t q and E^k based 
on this model yield more intelligible comparisons between 
data sets than those derived strictly by curve -fitting. 



TABLE 6.4.2 


BINDING ENERGY TEMPERATURES OBTAINED FROM At $ = T^exp^/k^) 

USING THEORETICAL x! BASED ON MOBILE AND LOCALIZED MODELS 

0 


r\ A t n r*rr - 

Mobile 

Localized 

DATA SET 

T o * sec 
{Theoretical ) 

E 0 /k , °K 
(Average) 

T 0 , sec 
(Theoretical) 

E o /k , *K 
(Average) 

Xe«*Cu 

2.5 x ID' 13 

2239 (± 3%) 

7. 6-9. 4 x 10" 17 

3155 (± 1%) 

Xe-»Ni A 

2.0 x 10" 13 

2294 (± 2%) 

3. 9-5. 2 x 10- 17 

3251 (± 0.3%) 

Xe-*-Ni A + D 

II 

2302 (±4%) 

3. 9-4.8 x 10 -17 

3239 (± 2 %) 

Xe ^Ni A + D + He + 

II 

2069 (±2%) 

3. 5-4. 4 x 10' 17 

2965 (± 0.4%) 

Xe-Ni A + D + Ar* 

II 

2080 (± 2 %) 

3. 5-4. 8 x 10" 17 

2996 (± 1%) 

Xe-*Ni A + D + Xe + 

II 

1971 (± 4%) 

3. 5-3. 9 x 10" 17 

2845 (± 1%) 

Xe^Ni B + D 

II 

2305 (± 2 %) 

3. 7-5. 2 x 10- 17 

3251 (± 1%) 

Xe-Ni B + D + Ar + 

II 

2132 (± 2 %) 

3. 3-4. 6 x 10" 17 

3026 (± 0.3%) 

Kr-*-Ni A + D 

1.9 x 10' 13 

1870 (± 1%) 

6. 6-7. 6 x 10“ 17 

2621 (± 1%) 

co 2 -*ni A + D 

2.0 x 10" 13 

2281 (± 3%) 

1.4-1. 7 x 10" 16 

3108 (± 1%) 


95% confidence intervals In average E /k based on n values for each 
data set are given in parentheses? 


oo 

VoJ 
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It can now be seen more clearly that small changes in 
binding energy produce order of magnitude changes in stay 
times. For Xe on Ni f the previously shown decrease in stay 
times by a factor of 10 to 100 due to ion bombardment is 
associated with a decrease in Ey^k of only about ?#. Simi- 
larly, the results for Kr on Ni indicate a binding energy 
some 20 % lower than that for Xe or C0 2 on Ni even though at 
comparable T g , the stay times might differ by a factor of 
1000, Also, a binding energy for C0 2 on Ni lower than that 

9 

for Xe and a higher for C0 2 than that for Xe are consis- 
tent with the qualitative behavior in Fig. 6,3.6 where the 
peak shift data were comparable in magnitude for the two 
gases. 


6.5 The Negative Dips in the Detector Signal 

The negative dip in the detector signal was discussed 
briefly in Sec. 6.1. For this dip to occur, the background 
component of the detector signal must decrease temporarily 
through some interaction with the beam pulse. It was stated 
earlier that this interaction most likely occurred between 
molecules adsorbed on the surface. To support this state- 
ment, the probability of a number of other possible 
interactions must be considered. First, it can be stated 
that the geometry of the apparatus was such that 
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interactions between beam and background molecules passing 
through the chopper slots and from there directly into the 
detector did not contribute to this dip. This was verified 
in separate experiments in which an additional aperture was 
placed in front of the detector such that no direct line of 
sight existed between the chopper slots and the detector 
ionization zone. Second, collisions between beam and back- 
ground molecules in the regions between the chopper, target 
and detector should be negligible. The mean free path of 
background molecules through gas at the density of the beam 
pulse is estimated to be at least 300 cm, much greater than 
the flight path lengths involved. Finally, the duration 
time between pulses was such that no anomalies due to pulse 
overlap should have occurred. For example, if the t 2 data 
are a valid upper limit for the mean stay time for this 
type of signal, it is estimated that no greater than 5% of 
the molecules from the preceding beam pulse remained on the 
surface at the beginning of a new pulse even at the largest 
t 2 encountered in these experiments. Thus, barring some 
unforeseen behavior, it seems that the beam-background 
interaction did occur most likely on the target surface. 

To estimate the effects of these dips on the data 
reduction, an empirical model was developed. It was 
observed that the position of the "minimum” occurred near 
the position of the expected maximum for the detector signal 
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for zero stay time. This suggested that the observed signal 
might be modeled by the following superposition of two 
signals i 


S(t,r,K) « S(t,r) - KS ( t , 0 ) (6.5.1) 

where K is an adjustable constant. The function 
represents the contribution to the detector signal from the 
desorbing molecular beam pulse. The term, -KS(t,0), repre- 
sents the decrease in the contribution to the detector 
signal from the background gas. For given t, temperature, 
and gas species, the functions S(t, < r) and S(t,0) were 
assumed to be given by Eq. (4.1.6). The combination of ~r 
and K was found which gave the best fit to both the measured 
peak shift and the ratio of the positive to negative detec- 
tor signal amplitudes S max /|S min |. Two typical comparisons 
of measured to modeled signals are shown in Fig, 6.5,1. In 
the first example, Fig, 6.5.1(a), the agreement is excellent 
for t >t_, but the model's minimum occurs at a later time. 

In Fig, 6.5.1(b), the predicted and measured minimum posi- 
tion are in slightly better agreement, but clearly the model 
cannot account for the observed exponential-like rise in the 
negative portion of the signal. Signals of the shape shown 
in Fig. 6.5.1(h) imply that the stay times are logg compared 
to the measurement times and therefore only the background 





S(t) # arbitrary units 
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dip is observable in the signal. The majority of the 

recorded signals showing the negative dip were of the general 

shape shown in Fig. 6.5.1(a). The appearance of this type 

of signal depended on the magnitude of the stay time. It 

-4 

did not appear for peak shifts less than 10 sec for any of 
the test gases. 

For one set of detector signals, values of t were 
deduced using Eq. (6.5.1) for those signals showing the dip. 
The result of these calculations is shown in Fig. 6.5.2. It 
is seen that -r deduced in this manner shows better agreement 
with Tg than the values given in Fig. 6,1.2. The new analy- 
sis seems still to overpredict the stay time at large t 
where 7*2 should be a reasonable upper limit for r. However, 
this overprediction is not as great as that shown previously 
in Sec. 6.2. 

Obviously, the empirical model gives little insight into 
the physical mechanisms behind the observed behavior. It 
seems that dips are indeed due to temporary reduction in the 
background-component desorption rate caused by the beam 
pulse. Such a reduction indicates that the effective bind- 
ing energy increases with coverage. The fact that the 
magnitude of the dip minimum was observed to increase with 
decreasing temperature (increasing coverage) also supports 
this conclusion. Such an increase in binding energy has 
been observed by Baker and Fox^-^ 


for xenon on nickel. 




Figure 6.5.2 Example of Stay Time Obtained from 
Peak Shift Based on Empirical Model 
for Negative Dip, Xe-*-Ni A 
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However, since the maximum coverage of beam gas species 
(beam pulse plus background component) is estimated to be a 
small fraction of a monolayer, it is difficult to see how 
interactions between beam-species atoms lead to the observed 
dip. Of course, there is a high probability that a signifi- 
cant coverage of other background species exists in this 
temperature. Suffice it to say that, in the absence of 
supporting information concerning the surface conditions, 
this behavior cannot be completely explained. 
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CHAPTER 7 

SUMMARY AND CONCLUSIONS 

It has been demonstrated that a molecular-beam time-of- 
f light approach is potentially powerful for the study of 
physical adsorption under non-equilibrium conditions. Only 
under non-equilibrium conditions can information about 
adsorption and desorption rates and probabilities be obtained. 
In particular, it has been shown that a mean stay time for 
molecules on the surface is a measurable quantity. It has 
also been observed that pronounced effects on the desorption 
pulse can occur due to what are apparently adsorbate- 
adsorbate interactions. 

A model for the detected time-of-flight signal behavior 
has been developed assuming a first-order desorption process. 
This model predicts that, for a mean stay time v comparable 
to or less than the mean transit times for the experiment, 
the detected signal is delayed by an amount approximately 
equal to this stay time. Therefore, the stay time can be 
estimated from the shift in the signal peak. For stay times 
much greater than these transit times, the model predicts 
that the signal decays as exp(-t/r) and r can be estimated 
from this decay rate. For intermediate values of t, no 
simple relationship between signal shape and t exists. The 
shutter open time t g has been shown to have negligible 
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effect on the peak shift due to stay time, but unless r»t g , 
it may have significant effect on the decay rate. 

From the experimental data, mean stay times have been 
estimated for the adsorption of Xe, Kr and C0 2 on copper 
and nickel surfaces as a function of surface temperature. 

The minimum detectable stay time for these experiments was 

_ c; _ c _ii 

found to be about 10 sec. For t between 10 J sec and 10 

sec, t was estimated from the measured peak shift. However, 

—ll 

for r >10 sec, interpretation of the peak shift in terms 
of the stay time becomes progressively less certain primari- 
ly because of a peculiar negative dip in the detector 
signals which is attributed to adsorbate-adsorbate inter- 
actions. 

The mean stay time has been found to be relatively 
independent of the surface materials and of the pre-vacuum 
surface preparations used in these experiments. This could 
be due to obscuring of the bulk material by a layer of 
physisorbed gases reaching the surface through the back- 
ground gas as well as chemisorbed gases present prior to the 
experiments. 

Ion bombardment of the nickel surface during Xe adsorp- 
tion has been found to reduce the mean stay time at fixed T . Xe 

s 

bombardment reduced the stay time typically by a factor of 
about 100 while Ar and He ion bombardment each reduced the 
stay time typically by a factor of 10. It is proposed that 
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this reduction results from burial of the neutralized ions 
near the surface. For argon ion bombardment, ion doses 
corresponding to about 0.1 monolayers were sufficient to 
stabilize the stay time at its reduced value. 

Experiments were performed for Xe on Ni at fixed sur- 
face temperatures over a range of incident beam flux by 
varying the molecular beam source pressure. At the highest 
source pressure employed, the measured stay time and calcu- 
lated beam flux indicate a surface coverage of about 10’^ 
monolayers. Decreasing the pressure to where the flux was 
halved had little effect on the stay time. However, reducing 
the pressure by a factor of 10 (flux by factor of 3) caused 
the stay time to decrease substantially. 

The stay time for CO^ on nickel was found to be compar- 
able to that for Xe at the same surface temperature. Stay 
times of about the same magnitude were found for Kr on 
nickel only at temperatures about 20 % lower. 

The measured variation of the peak shift dt with sur- 

s 

face temperature was used to estimate adsorption binding 
energy and preexponential factors by fitting these data to a 
Frenkel-type relation, dt g = r' 0 exp(E 0 /kT s ). However, 
systematic uncertainty and scatter in the peak-shift data 
are such that values of *r and determined by direct 
curve-fit to this relation show large uncertainties. 

f 

Theoretical estimates of <r 0 based on a localized adsorption 
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model were found to be in reasonable agreement with the 

t 

values found by fitting. These values of *r are quite low 
sec) and are not representative of typical vibra- 
tional periods for adsorbed atoms. Because of the high 
uncertainty in the fitted -r and E c /k values, the localized 

I 

adsorption model estimates for 'r Q were used to obtain more 
reliable estimates of the relative values of E Q /k from the 
peak-shift data. The E Q /k value for Xe on nickel was thus 
found to be about 3250°K with a reduction to about 2850°K 
upon ion bombardment with Xe. For C0 2 and Kr on nickel, the 
values of E Q /k were 3100°K and 2600°K, respectively. 

Comparisons of measured signal shapes to those predicted 
by a first-order desorption model show reasonably good agree- 
ment at small stay times (low coverage). The measured 
signal peak shifts at small stay times behaved with surface 
temperature as expected. Significant disagreement between 
measured and predicted decay rates at long detector times 
had little effect on the peak shift behavior and were attri- 
butable to shutter function approximations in the model. At 
the longer stay times encountered in these experiments, the 
adsorbate -adsorbate interactions, proposed to account for 
signals exhibiting the negative dip, obviously violate the 
assumption of first-order desorption. Still, at long mea- 
surement times in the detector signal, the signal decay rate 
was indicative of first-order desorption. 
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By far the most serious shortcoming of the present 

experiments is the uncertainty with regard to the nature of 

the surfaces. The relatively high partial pressures 
-8 

( w 2 x 10 torr) in the test chamber of species other than 
that in the molecular beam could have led to significant 
adsorption of these species. The presence of such adsor- 
bates may well be the cause for the unexpected negative dips 
seen in the detected desorption signals. Uncertainty in the 
constancy and nature of the surface conditions also makes 
the interpretation of the effects of parameter variations 
correspondingly uncertain. 

To improve on these experiments, the following recom- 
mendations are offered t 

(1) employ an ultra-high-vacuum test chamber 

capable of attaining residual pressures of 
-1 0 

about 10 torr (with sufficient pumping 
capacity to keep the pressure below about 
10“ 6 torr with the molecular beam on)i 

(2) use single crystal targets » 

(3) use in-situ cleaning techniques that allow 
the production of atomically clean surfaces 
(if ion bombardment is used, have the capa- 
bility of heating the target sufficiently to 
remove buried ions and to anneal the surface )» 
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(4) employ in-si tu analysis of surface, such as 
Auger spectroscopyi 

(5) maintain ability to vary incident flux in 
order to study the effect of adsorbate- 
adsorbate interactions; 

(6) investigate the effect of ion burial on 
adsorption (if the effect is as pronounced as 
it appears, it is potentially useful for 
altering surface adsorptive characteristics)! 

(7) determine the physical mechanism which causes 
the negative signals (this might provide an 
approach to studying adsorbate-adsorbate 
interactions) . 
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APPENDIX I 

PREDICTED STAY TIMES FOR ADSORBED GASES 


It was shown in Chapter 2 that the stay time is related 

to the number ratio- N /N , bv 

a g 



/ 2jra 

k* v / N «\ 


T - 

2 

— — I 



\ kT ) 

A_ \N / 
8 g 

(AI.l) 


for a steady-state adsorption process with unity adsorption 
probability. Stay time relations will now be derived for 
non-interacting;, indistinguishable molecules in equilibrium 
with a surface for some specific models of the gas-surface 
interaction. 

Mobile harmonic-oscillator model . In this model, the 
adsorbed gas is assumed to move freely in the two dimensions 
parallel to the surface. Molecules are assumed to be bound 
to the surface with energy -E 0 (the zero-energy level is 
taken to be at the beginning of the continuum state) and to 
be vibrating perpendicular to the surface as simple harmonic 
oscillators. At equilibrium, for non-interacting, indistin- 
guishable molecules, the ratio, N /N , is given by 

^ g 
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where Q a and are the single-particle partition functions 
for the adsorbed phase and gas phase, respectively (for 
example, see Reif^ 3 ^). If the gas is of sufficiently low 
density, the classical limit 

= V (m g kT/27Tfi 2 ) 3//2 (AI.3) 

may be used. The partition function Q may be written as 

ci 

Q a = Q a rQ a ibexp(E o/ kT) (AI.4) 

where Q a = translational partition function* 
vib 

Q a = vibrational partition function. 

The classical limit 

«a r = “gM/Sa* 2 ( A i. 5 ) 

may be used for the translational part. For the vibrational 
part, however, it is appropriate to use the quantum- 
mechanical result for a one-dimensional, harmonic 
oscillator. 



(AI.6) 
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where u) Q = frequency of vibration. (For evaluation of 
partition.. functions for adsorbed gases, see Reif^-^ or 
Knuth^-^.) Combining Eqs. (AI.1)-(AI.6), the stay time is 
found to be 

ficd> /2kT 

T = — o 

sinh(fio> 0 /2kT) 

where v Q = 2 n/w 0 is 'he vibrational period. 

Localized, harmonic oscillator model . In this case, 
consider the molecules to be bound so strongly in localized 
sites that their translational freedom is lost entirely. 

By accounting for the number of ways N atoms can be dis- 

Ck 

tributed over N m sites (N m corresponding to monolayer 
coverage), it can be shown that 


T 0 exp(E Q /kT) (AI.7) 



(AI.8) 


For the localized case, Q a =1. Treating the molecules as 
three-dimensional harmonic oscillators with frequencies cj ± 
and for vibrations perpendicular and parallel to the 
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surface, respectively, the vibrational partition function is 

(^ 7 ) 

given by Wf ' 



The result for the stay time then becomes 


(AI.9) 




exp(E 0 /kT) 

m (kT) 2 ' 

c> 


2sinh{ 

^2kT/ 
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2sinh| 

ffiwA 
{ 2kT / 


(AI.10) 


Calculation of the vibrational frequencies . Character- 
istic frequencies for the surface atoms can be estimated 
from the Debye temperature of the solid, i.e., 

<*> = kflyti (AI.ll) 

For surface atoms, a better estimate is obtained if a "sur- 
face Debye temperature" such as that defined by Goodman 
is used, viz., 

^D, surf = % (AI.12) 


where p is the ratio of mean square displacement of a 
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surface atom to that of a bulk atom. Values of p have been 
calculated by Allen and deWette^^ and p ~ 2 is typically 
found. (There exists considerable variation in theoretical 
values of p as found in the literature. The quoted value 
represents the author's best estimate for both inert gases 
and metals based on the literature.) Assuming - 

= kSjy^h, Eqs. (AI.,7) and (AI.10) become 


for mobile case » 


T 


Op/k T 

sinhffljjAT) 


T o exp(E Q /kT) 


(AI.13) 


for localized case i 

4mh 2 <7 ra (ej/2T) 2 
m^k# D ^sinhC^o/^jJ^ 


T 0 exp(E 0 /kT) 


(AI.14) 


where a = N /A . Both relations can be written in the 
m nr s 

general form 

T = Vexp(E 0 /kT) (AI.15) 

These results differ from that of Frenkel only in terms of 

« 

the preexponential factor t q . 

Sample calculations for Xe on Ni . Using = 413 °K 
for nickel ^ 28 ^ and a m = 4.8 x 10 lif atoms/cm 2 for the 
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monolayer coverage of Xe (29) p the factor t q was calculated 
at T = 100 °K. The results arei 

T 0 ~ 2.3 * 10 * 3 sec (Frenkel result )» 

T o, mobile ~ ^ sec * 

^o, local “ 3,6 x 1°" 17 sec. 

It is seen that the pre exponential factor for the localized 
case may be orders of magnitude below that for a mobile 
adsorbed layer. The factor t q is associated with the loss 
of entropy upon adsorption*^ and thus, in the localized 
case, loses the meaning of a vibrational period. 
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CALIBRATION OF THE DETECTION SYSTEM AND CHARACTERIZATION 
OF THE MOLECULAR BEAM SOURCE 

Ionization Detector . The flow-through-type ionization 
detector used in this study is sensitive to molecular density 
in the ionization region. The geometry of the system is 
similar to that used by Hagena, Scott, and Varma^ 2 ^, It is 
different in that a 900 Oe magnetic field is used to confine 
the electron beam (compared to about 350 Oe in Ref. 26) and 
in that 1 % thoriated (rather than pure) tungsten filaments 
are used, allowing operation at lower filament temperatures. 
The performance of the detector was determined in a small 
bell- jar vacuum system. The detector ion current was meas- 
ured as a function of chamber pressure (indicated by a 
commercial, calibrated, hot-cathode ionization gauge), 
filament potential, anode potential and ion collector 
potential. A plot of measured ion current versus chamber 
pressure is shown in Fig. AII.l. The potentials listed in 
this figure are approximately those which gave maximum 
sensitivity for a given electron emission current. The ion 
current should be linearly proportional to the gas density 
within the detector’s active zone. The observed variation 
with pressure is approximately linear with the standard 
deviation of all individual measurements from a least- 
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Chamber Static Pressure t torr 


Figure AXI.l Variation of Detector Ion Current 
with Chamber Pressure, N_ Gas, 
Emission Current = 10 mA. 
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squares straight line being about 9 %. This deviation could 
well be due in part to lack of precision in the instrument 
used for measuring chamber pressure. 

For the stay —time experiments* the signal is a time- 
varying pulse and the detector titne response must be known. 

A straightforward way to evaluate this response is to use a 
pulsed molecular beam for which the detector density-versus- 
time is accurately known. For An orifice operating at a 
sufficiently large Knudsen number Kn (the ratio of mean free 
path in the source to orifice diameter ) and pulsed with high 
time -of -flight resolution^ 2 ^, the signal for a density 
sensitive detector is of the form 

S(t) « ^ exp(-0 o 2 L 2 /t 2 ) (All. 1 ) 

where t * flight time, 

0 o = vW2kT 0 , 

L = flight path length. 

Detector signals were obtained using a 0.04 cm diameter 
single-orifice source with nitrogen gas. The results are 
shown in Fig. All. 2 for two values of Kn. For comparison, 
the corresponding Maxwellian distribution from Eq. (AII.l) 
is plotted in this figure. The agreement between measured 
and predicted results is very good for Kn - 5.7 but not for 




(a) Kn = 5.7 


Figure All. 2 


Comparison of Measured to Predicted Incident 
Beam Time-of-Flight Distributions for a 
Single-Orifice Source, N 0 Gas 



(b) Kn = 0.6 


Figure All. 2 concluded 
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Kn=0.6. At Kn = 5.7, the measured most probable time t m 
was within t 1% of that predicted by free-molecular theory* 
the measured signal time-width at %S(± m ) for Kn = 5.7 was 
within i 6% of the theoretical value. The assumption of 
free-molecular effusion through the orifice at Kn = 0.6 is 
no longer justified and the actual disagreement is similar 
to that generally observed in this Kn range. 

As a result of a series of measurements of this type 
and their comparison to theory and previous experiments, it 
can be stated that the detector signals were representative 
of the time variation of density within the detector's 
ionization region. 

Multichannel Molecular Beam . Following the tests 
described above, the detector system was used to characterize 
the molecular beam produced by the multichannel source. For 
the stay-time experiments, the properties of interest are the 
intensity and the speed distribution on the beam axis. Also, 
since the small diameter channels of the multichannel source 
were highly susceptible to blockage by dust particles, it 
was important to periodically determine the intensity during 
the stay-time experiments to ensure that no blockage 
occurred. Since direct measurements of the intensity 
required removal of the target and relocation of the detec- 
tor, the centerline intensity I was estimated indirectly by 

I 0 = X2/tt (All. 2) 



114 


where X is the peaking factor defined by Jones, Olander, and 
(k\ ) 

Kruger' ' and l is the total molecular flow rate. By 
definition, the peaking factor is the ratio of the center- 
line intensity to that from an effusive source emitting with 
a cosine distribution at the same total flow rate. This 
ratio was determined from intensity measurements using the 
effusive source described earlier and the multichannel source 
and was found to be about 6.6 for nitrogen. Although no 
measurements of x were made for other gases, x was estimated 
to be about 5 for Xe using Eq. (10) of Ref. 4l . 

The total flow rate through the multichannel source is 
related to the upstream number density n Q by the expression 

2 = n Q C (All. 3) 

where C is the source conductance. Knowing the source cham- 
ber pressure, temperature and volume, this conductance was 
determined by measuring the decay of the source chamber 
pressure with time with the gas supply shut off. The results 
of a typical measurement of the pressure decay is shown in 
Pig. All. 3. For xenon at a temperature of 295 to K, the con- 
ductance was found to be 8 x 10*^ i/ sec, within \% of that 
predicted by kinetic theory. This result was confirmed by 
several separate measurements made periodically during the 
course of the stay-time experiments. Using these results in 



Figure All. 3 Typical Measured Source Pressure Decay with Time 
Used for Multi channel -Source Conductance 
Determination, Xe Gas 
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Sqs. (All. 2) and (All. 3), the centerline intensity was cal- 
culated to be 4.5 x 10*^ atoms . 3 ter”* . sec”* for xenon. 

The corresponding centerline flux at a target 3.6 cm from 
the source exit is then about 3.4 x 10*-’ atoms . cm” 2 . sec”*. 
The uncertainty in the flux determined in this way is pri- 
marily due to the uncertainty in the peaking factor. The 
quoted value for Xe is believed to be accurate only to with- 
in a factor of 2. Whatever the factor, the flux was constant 
within ± 1 % for the experiments (except when source pressure 
was intentionally varied). 

Characterizations of N 2 and Ar beams produced by the 
multichannel source have been reported in a separate study 
A typical measured time-of-flight signal for N 2 taken from 
this study is shown in Pig. All. 4 where it is compared to 
the Maxwellian time-of-flight distribution predicted by 
Eq. (All. 1 ) . The measured distribution is narrower than the 
Maxwellian and results in a mean molecular energy some 
25 to 30 % higher than that for the Maxwellian distribution. 
This is attributed in Ref. 24 primarily to collisions at the 
exit of the multichannel array. Although no direct beam 
time-of-flight measurements were made for the gases used in 
the present stay-time experiments, these gases are not 
expected to produce significantly different distributions 
since the Kn for the present study is not much different 
from that for the distribution in Pig. All. 4. 




Figure All, 4 Comparison of Measured Incident Beam Time-of-Flight 
Distribution for Multichannel Source to Maxwellian 
Prediction, p = 10 torr, N~ Gas 



APPENDIX III 


DERIVATION OF THE EQUATIONS FOR THE SHAPE OF THE 
INCIDENT AND DESORPTION PULSES 

An approximation for the shape of the incident flux 
pulse F^t) is developed in this appendix. This approxima- 
tion is then used to predict the form of the desorption 
pulse for a first-order rate process. 

Incident Pulse . Consider a molecular beam with uniform 
flux F q such that the rate at which molecules with speeds 

between v and v + dv strike an element of area dA_ at a 

s 

particular downstream location is given by 

dF 0 = Cvf o (v)dv (AIII.l) 

In this expression, f“ 0 (v) is the incident beam speed distrib- 
ution and C is a constant of proportionality. Since only 
the timewise shape of the pulse is of interest, all factors 
affecting only the pulse amplitude will be included in the 
constant C. A shutter is now introduced between the source 
and target at a distance x^ from the target with a function 
g(t' ) representing the instantaneous fraction of molecules 
that pass through this shutter between t* and t’ + dt' 
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headed for dA,,. The flux incident upon dA_ at time t is then 
s s 





oo 


dv 

d(t-t* ) 


dt- 


(AIII.2) 


Substituting v = x^At-t' ) gives the result 

Fi<t) = F ° / 7 ^ 7 ? f °&) g(t ' ,dt ’ (Am - 3) 

- 00 

It is now assumed that the incident beam has a 
Maxwellian speed distribution, i.e., that 


f o(r) « t^expt-^Vt 2 ) (AIII.4) 

(Although the beam produced by the multichannel source does 
not have this distribution exactly, it will be shown that, 
for small x x , the shape of ? i (t) is only weakly dependent on 
f^Xj/t). ) 

The shutter function for a uniform beam of rectangular 
cross section with width equal to the chopper slot width is 
the triangular function 

g(t) * 1 + t/t s for -t s <t<0 

= 1 - t/t g for 0 < t < t g (AIII.5) 

= 0 otherwise. 



120 


where 2t g is the basewidth of the triangle. It can be shown 

that, if t g is allowed to become very large compared to any 

transit times and stay time, the resulting detector signal 

will approach the shape of g(t). This fact was used 

experimentally to measure the effective value of t (see 

s 

Sec. 5.2). In addition, based on these measurements, the 

shape of the true g(t) was approximately determined and 

found to be nearly Gaussian, i.e., g(t) oc exp(-Bt^/t ^), 

s 

Thus, for the present experiments, Eq. (AIII.5) describes 

the shutter function only approximately, (It can further 

be shown that the detector signals predicted in Chapter 4 

are sensitive mainly to the effective width and not to the 

exact shape of g(t), even for t w 0 x„. ) 

S S 2 

Substituting Eqs. (AIII.4) and (AIII.5) into Eq. (AIII.3) 
and integrating yields 

for -t < t < Of 

s 

O s 

- iV5Wt„)[i - «*(*&-)] 

s 
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for 0 < t < t _» 

s 

= i(l +t A s )exp[-(^-) ] - (t/t s )exp[-(<J/t) 2 ] 

+ i^/t g )[erf(^| i -) - 2erf (i/t) + l] 

s 

2 

i(l+t/t s )exp[-(^-) ] - (t/t g )exp[-(d/t) 2 ] 

f s 

- *<l-t/t s )exp[-(^-) 2 ] 

s 

+ iV5WV[erf(^) - 2erf (i/t) + erf(^)] 

(AIII.6) 

where 6 - Eq. (AIII.6) was evaluated for several 

values of t g /d covering the range of values expected in the 

stay-time experiments. The result for t /j = 10 is shown as 

s 

the solid curve in Fig. AIII.l. The fact that the calculated 
shape very nearly corresponds to the assumed triangular 
shutter function shows that the incident beam speed distrib- 
ution is not of much importance for values of t /i as large 
as this. To first approximation, only the mean transit time 


for t g < t* 
F ± ( t ) 
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is important. The approximation 


= i + <t-t 01 )/t a 
" 1 - < t - t ol ) / t s 


for -t g < (t-t Ql )< 0 
for 0 < (t-t Ql ) < t B 
otherwise 

(AIII.7) 


where t Ql = 6 /\pZ, is also shown in this figure for 
comparison. Since t /$ was always greater than 10 in these 

o 

experiments, the approximate expression for F i (t) is used in 
interpreting the recorded detector signals. 

Desorption Pulse . For first-order desorption, the 
adsorbate coverage is related to F i (t) through the following 
integral expression i 


(t) = rF exp(-t/r) 


(T 


) 


-exp(t'/r)dt' 


-00 


(AIII.8) 


Substituting Eq. (AIII.7) into the above and integrating 



a 24 


gives 


for -t g < (t-t Ql )< °» 


<y (t) f /t+t -t \ 

irP = 1 + (t - t ol" r)/t s + <' r /V exI> -( T ) 

o L ' 


for 0 < (t-t Ql ) < t g j 
<r (t) 

tF = 1 - “ 2(r/t g )exp 


-PM 


+ (r/t g )exp 




for 0 < (t-t Ql ) » 


<7 a ( t ) 

= (v/t 3 )exp 


+ (r/t g )exp 


[-m] • 
[-pmm 


2(r/t g )exp 


-phi 


(AIII.9) 


The quantity, 'r? 0 , is simply the surface coverage that would 
be obtained if the beam was not pulsed (i.e., chopper removed). 
(It should be noted that, if the actual beam cross-section 
is greater than the chopper slot-width, the maximum coverage 
for the pulsed case will differ from tF q by some fixed 
geometrical factor. ) 
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The ratio, <? a (t )/ tF o , is plotted versus the non- 
dimensional time t/d in Fig, AIII.2 for several values of 
r/t s ' * ?or = 0,1, the shape of the desorption pulse 

(desorption flux equal to a & ( t)/ r ) is essentially the shape 
of the incident pulse only delayed by an amount -r/d. On the 
other extreme, for t/ t g = 5, the desorption pulse decays at 
large t/d approximately as exp(-t/V). 

An important feature of the predicted desorption pulse 
is shown by examining the behavior of a (t) at large x/t 
It can be shown that 


lim 

T/t-O, 


°a (t) 

TP o 


at 

t=t 


m 


T 


(AIII.10) 


Thus for T»t g , the maximum surface coverage is linearly 
proportional to the effective shutter open time. 











APPENDIX IV 

TABULATION OF STAY-TIME PARAMETERS 

The stay-time parameters obtained from the recorded 
detector signals for eaph of the data sets are Listed in 

Table AIV.l. Although the values of r calculated from At 

-4 S 

are in error for 4t^ > 10 sec, they are listed here for 

reference and to indicate the relation of ~r to the measured 

peak shift. 
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TABLE AIV.l 


TABULATION OF STAY-TIME PARAMETERS 


DATA SET 


Xe - Cu 


Xe - NiA 


Xe * NiA+D 


Xe - NiA+0+He + 


Xe - NiA+D+Ar + 


Xe - Ni A+D+Xe + 


Xe + NiB+D 


Xe ^ NiB+D+Ar + 


Kr -*• NiA+D 


CO^ -► NiAfD 


T $ , °K 

At s , usee 

usee 

t, usee 

116,7 

25 

350 

25 

114.1 

58 

391 

60 

110,4 

151 

517 

234 

108.1 

308 

706 

996 

105.3 

738 

1548 

9704 

102.4 

1450 

3774 

> 10 4 

122.7 

14 

408 

14 

119.7 

32 

431 

31 

117.1 

55 

474 

56 

113.7 

no 

501 

140 

111.4 

207 

822 

426 

108.1 

538 

1637 

4069 

104.9 

1275 

4289 

> 10 4 

117.0 

27 

367 

27 

114.1 

74 

415 

80 

111.1 

218 

970 

471 

108.2 

730 

1003 

9624 

105.1 

1335 

-* 

> 1 G 4 

111.2 

16 

382 

16 

108.3 

32 

380 

32 

105.3 

66 

430 

71 

102.3 

139 

635 

203 

99.1 

404 

1115 

1811 

117.0 

11 

359 

11 

114.1 

12 

354 

12 

111.2 

19 

383 

23 

108.3 

31 

451 

31 

105.2 

77 

407 

85 

102.4 

176 

618 

306 

98.6 

602 

1094 

5122 

105.0 

20 

- 378 

20 

102.5 

46 

426 

47 

98.7 

126 

636 

173 

122.6 

15 

688 

14 

119.8 

30 

730 

30 

117.0 

45 

817 

45 

114.2 

119 

751 

157 

111.0 

245 

1180 

602 

108.3 

594 

2040 

5348 

105.2 

1404 

2611 

> 10 4 

113.9 

18 

370 

18 

111.1 

29 

344 

29 

108.1 

62 

400 

65 

104.9 

123 

552 

166 

101.8 

288 

946 

830 

98.7 

736 

1786 

9077 

95.4 

1738 

5116 

> 10 4 

98.6 

20 

304 

21 

97.3 

37 

314 

37 

95.6 

55 

343 

58 

93.8 

87 

409 

102 

91 . 9 

164 

737 

292 

122.5 

15 

190 

14 

119.8 

24 

201 

23 

117.0 

48 

235 

48 

114.0 

1-12 

373 

152 

111.1 

246 

617 

983 

108.1 

665 

1273 

> 10 4 



